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摘要 ”作为 聚变 反应 堆 候 选材 料 ， 低 活化 铁 素 体 钢 〈 马 氏 体 钢 ) 的 主要 成 分 是 Fe 和 Cr。 在 实际 工 况 下 ， 


z 
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氨 会 进入 材料 中 ,产生 微观 缺陷 ， 从 而 引起 辐 照 损伤 ， 影 响 材料 的 抗 辐 照 能 力 〈 所 脆 现 象 ) 


聚变 反应 ] 
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及 其 它 物理 性 能 。 本 文通 过 第 一 性 原理 计算 ,模拟 得 出 反应 堆 主 体 材料 Fe-Cr 合金 的 界面 ， 在 不 


的 蔡 代 


位 和 间隙 位 迭 杂 He 原子 ， 对 其 结构 优化 后 计算 每 种 结构 的 形成 能 ， 得 出 He 原子 在 Fe-Cr 界面 不 同位 置 处 


SS 


的 稳定 性 。 与 此 同时 , 通过 对 体积 的 计算 得 出 He 在 靠近 界面 的 Fe 八 面体 间隙 位 时 会 引起 较 大 的 晶 格 畸变 ， 


根据 电子 态 密度 分 析 ， 推 断 是 He. Fe. Cr 三 种 原子 在 特定 位 置 相互 杂 化 影响 的 结果 ， 为 下 一 步 研究 核 级 


不 锈 钢 微 观 缺陷 及 其 物理 化 学 性 质 提 供 理论 依据 。 
关键 词 ”Fe-Cr 合金 ;第 一 性 原理 ， 氨 ; 界面 
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Abstract: As a candidate material for fusion reactors, the main components of low activation ferrite steel 
(martensitic steel) are Fe and Cr. In actual working conditions, helium produced by fusion reactions will enter the 
material, causing microscopic defects and radiation damage, affecting the material's resistance to radiation (helium 
embrittlement phenomenon) and other physical properties. This article simulates the interface of Fe Cr alloy, the 
main material of the reactor, through first principles calculations. He atoms are doped at different substitution and 
interstitial sites, and the formation energy of each structure is calculated after optimizing its structure. The stability 
of He atoms at different positions on the Fe Cr interface is obtained. At the same time, by calculating the volume, 
it was found that He would cause significant lattice distortion at the Fe octahedral gap near the interface. Based on 
the analysis of electronic density of states, it is inferred that this is the result of the hybridization of He, Fe, and Cr 
atoms at specific positions, provide a theoretical basis for the next research on the microscopic defects and 


physical and chemical properties of nuclear grade stainless steel. 
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根据 我 国 《中 国 绿色 能 源 发 展 纲 要 〈2023-2033)》 和 “《“ 十 四 五 ”能 源 规划 》， 核 能 正 逐 
渐 成 为 我 国 重要 的 绿色 能 源 ， 要 积极 有 序 发 展 核能 ， 而 核反应 堆 材 料 对 于 氢 所 的 抗 辐 照 能 力 
其 一 直 是 制约 和 限 值 核 能 发 展 的 因素 之 一 。 作 为 未 来 核 聚变 反应 堆 候选 材料 中 的 热门 材料 一 
一 低 活化 铁 素 体 / 马 氏 体 钢 的 主要 成 分 就 是 Fe-9Cr 合金 "， 其 国内 外 主要 成 分 如 表 1 所 示 。 
而 在 工 况 下 ， 高 能 粒子 辐 照 会 让 材料 内 部 产生 缺陷 , 包括 空位 及 相关 原子 ， 其 中 主要 包括 空 
位 "、 间 隙 、 位 错 及 位 错 环 等 复杂 结构 ”"， 这 些 复 杂 结 构 会 破坏 材料 的 内 部 构造 ， 从 而 导致 
影响 材料 的 抗 辐 照 、 抗 蠕 变 等 能 力 。 因 此 , 本 文选 择 使 用 核反应 堆 中 最 常见 的 He JR FE AB 
杂 原 子 ， 模 拟 其 在 Fe-Cr 合金 界面 时 的 情况 ， 从 而 进行 计算 分 析 。 

表 1 国内 外 低 活化 铁 素 体 / 马 氏 体 钢 的 主要 成 分 [4-9] 


Tab. 1 Main components of low activated ferrite / martensitic steel all around the world'*”! 


国家 型 号 主要 合金 元 素 微量 合金 元 素 
欧洲 EUROFER97 Fe (基体 ) ,9Cr W, V, Ta 

美国 9Cr2WVTa Fe (基体 ) ,9Cr 2W，V，Ta 

日 本 F82H Fe (基体 ) ,7.5Cr 2W, V, Ta 

日 美 联 合 JLF-1 Fe (4%) ,9Cr 2W, V, Ta 

日 美 联 合 JLF-2 Fe (4%) ,9Cr 2W, V, Ta 

日 美 联 合 JLF-3 Fe (基体 ) ,7Cr 2W，V，Ta 

中 国 CLAM Fe (基体 ) ,9Cr 1.5W, V, Ta 

中 国 CLF-1 Fe (基体 ) ,8Cr 1.5W, M, V, Ta 
中 国 SCRAM Fe (26/6) ,9Cr 2W, V, Ni, Ti 


核 材 料 领 域 中 ， 研 究 人 员 针 对 于 Fe-Cr 合金 ,主要 
合 


没 等 方法 展开 了 微观 缺陷 研究 。 有 研究 表明 ，Fe-Cr 


过 扫描 电镜 、 透 射电 镜 、 正 电子 渣 
中 


Cr 原子 的 适量 添加 会 减弱 Fe-Cr 合金 中 位 错 及 空洞 的 形成 ””。 对 Fe-9%Cr 合金 进行 退火 实 
验 ， 当 温度 上 升 至 500°C, 较 大 的 微观 缺陷 已 经 逐渐 消失 ， 当 温度 继续 升 高 直至 800°C 后 ， 
Cr 原子 开始 聚集 并 形成 Cr 析出 相 ， 最 终 溶解 ”。Fe-Cr 合金 作为 核反应 堆 第 一 壁 的 基体 材 
料 , 研究 人 员 通 过 离子 注入 机 ， 电 化 学 等 多 种 方法 ,针对 其 微观 缺陷 内 部 的 氨 迁 移行 为 开展 
了 大 量 的 实验 研究 “”， 对 于 Fe-Cr 合金 中 的 氨 行 为 推断 入 总 结 如 下 : 单个 或 多 个 氮 原 子 进 
入 材料 后 容易 与 空位 和 其 他 氨 原 子 相 结 合 ， 从 而 聚集 形成 氨 泡 , 这 一 现象 在 实验 中 能 很 好 的 
体现 出 来 ,同时 在 辐 照 的 作用 下 会 伴随 着 Frenkel 缺陷 对 产生 。 当 越 来 越 多 的 氨 原 子 进 入 后 
会 引起 品格 畸变 ,最 终 与 空位 不 断 吸引 形成 团 徐 ， 最 终 产 生肉 眼 可 见 的 氨 泡 ， 影响 材料 的 物 
理化 学 性 质 。 


本 次 研究 通过 采用 密度 泛 函 理论 ， 通 过 第 一 性 原理 计算 ， 可 以 高 精度 的 计算 得 出 He 原 


子 在 Fe-Cr 界面 不 同位 置 的 稳定 性 ， 从 理论 上 得 出 He 


原子 在 不 同位 置 时 对 Fe-Cr 界面 品格 


参数 的 影响 ， 为 核 级 不 锈 钢 界面 处 氨 迁 移行 为 的 研究 提供 理论 基础 。 


1. 计算 方法 


第 一 性 原理 又 名 从 头 计算 ,计算 特点 是 只 考虑 粒子 之 间 的 作用 原理 ,通过 量子 力学 的 方 
法 去 求解 薛 定 齐 方程 ， 从 而 计算 得 到 材料 物理 化 学 性 质 的 方法 ,计算 过 程 中 不 含 任何 经 验 参 


数 。 本 次 使 用 第 一 性 原理 计算 软件 是 维也纳 从 头 计算 模拟 包 VASP) ”。 


2. 结构 优化 


本 次 优化 的 结构 模型 设置 为 Fe-Cr 界 面 ， 属 于 典型 


的 体 心 立方 结构 。 对 于 氨 原 子 挫 杂 的 


位 置 设 置 为 蔡 代 位 和 间隙 位 ,其 中 间隙 位 又 分 为 八 面体 间隙 和 四 面体 间隙 。 本 论文 使 用 VASP 
软件 对 Fe-Cr 界 面 进行 结构 优化 ， 优 化 方法 采用 GGA 交换 关联 泛 函 〉 和 PAW (级 加 平面 波 ) 


25。 另外 设置 TSIF=3，ENCUT=400 eV， 直 到 优化 至 每 


个 原子 受 力 二 0. 001 eV。 优 化 结束 后 


继续 进行 静态 计算 ，IBRION=-1， 得 出 界面 处 原子 分 布 模型 (图 1)。 


1 Fe-Cr 合金 界面 的 原子 模型 (1,2,3,4 分 别 代表 平行 于 xoy F 


面 的 切 Hl, Fe 原子 用 5 TERI, Cr 原子 用 


st 


蓝 色 表示 ) 


Fig. 1 Atomic model of Fe-Cr alloy interface (1, 2, 3, 4 represent the sections parallel to the xoy plane, respectively, 


yellow represents the Fe atom and blue represents the Cr atom. ) 


Ale, Cr 原子 用 蓝 色 表 示 ，Fe 原子 用 黄色 表示 ， 切 线 1、2、3、4 分 别 表 示 He 原子 


所 在 的 平面 。 另 外 ，T (Tetrahedron) 表示 He 原子 位 于 四 面体 间隙 位 ，0 (Octahedron) 


表示 He 原子 位 于 八 面体 间隙 位 、S (Supersede) 表示 He 原子 位 于 替代 位 ， 为 了 进行 对 
比 ， 另 加 入 空位 : V Vacancy) 表示 空位 在 替代 位 《图 中 某 个 Fe 或 Cr 原子 消失 ， 其 原 
所 在 的 位 置 没 有 任何 原子 )。 比 如 : T3 就 表示 He 原子 位 于 3 号 切面 上 ， 并 处 于 由 一 个 Cr 
原子 和 三 个 Fe 原子 所 构成 的 四 面体 的 中 心 。 在 结构 优化 的 基础 上 ， 分 别 把 He 原子 添加 
到 以 上 位 置 处 ， 再 对 其 进行 结构 优化 ， 优 化 参数 与 之 前 相同 ， 优 化 后 再 进行 一 次 静态 计 
算 ， 并 利用 以 下 公式 计算 其 形成 能 ”: 


EFormation = Ere. Cr. He 一 Exe. Cr 一 EHe (1) 


E Formation = Ere. Cr. He 一 Ere. Cr EHe + Esingle Fe/Cr (2) 


FHA, En o ne ANBAR A He 原子 的 Fe-Cr 界面 体系 的 总 能 量 ，E 6 表示 Fe-Cr 界 
HKR, En 表示 单个 He 原子 在 一 个 相对 无 穷 大 空间 内 的 ， 

示 被 替代 的 单个 Fe 或 Cr 原子 的 静止 能 量 。 在 经 过 相关 计算 过 后 ， 得 出 来 的 He 在 Fe-Cr 
界面 各 个 位 置 的 形成 能 如 图 2 所 示 。 
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图 2 He 和 空位 在 Fe-Cr 合金 界面 不 同位 置 处 的 形成 能 


Fig. 2 Formation Energy of He and Vacancies at Diffe rent Positions of Fe-Cr Alloy Interface 


3. 形成 能 


MBI VI, V2, V3 及 V4 的 形成 能 可 以 得 


UK 


Fe-Cr 界面 中 少 了 一 个 原子 ， 导 致 


co 


单 空 位 产生 时 ， 其 结构 是 不 稳定 的 。 相 比 于 Cr 原子 处 出 现 空位 ， 空 位 


性 也 相对 最 高 )。 换 句 训 


6 来 讲 ， 当 Fe-Cr 界面 中 出 现 单个 空位 时 ， 空 位 


更 倾向 于 在 远离 界 


AA) Fe 原子 位 置 形成 V4 位 置 表示 远离 界面 的 Fe 原子 位 置 ， 其 形成 能 相对 最 小 ， 稳 定 


更 倾向 于 在 远离 界 


根据 图 中 He 原子 


的 形成 能 可 以 得 出 ， 氨 原子 在 四 面体 间 际 位 时 相对 最 稳定 ， 其 次 是 在 


蔡 代 位 时 相对 最 不 稳定 。 在 这 之 中 ， 氨 原子 在 Tl 位 置 处 的 形成 能 最 低 ， 在 S1 位 置 处 的 


形成 能 最 高 。 由 此 可 以 得 出 ， 当 单一 He 原子 注入 Fe-Cr 界面 后 ,会 首 


MAY Cr 四 面体 间 阶 位 中 ， 如 有 多 个 He 原子 进入 ， 会 进入 其 他 四 面体 
界面 的 Fe 四 面体 间隙 位 都 填 满 后 ， 才 会 进入 八 面 体 间隙 位 中 ,一 直到 


MAY Fe 原子 位 置 处 形成 ， 即 四 周 都 是 Fe 原子 的 位 置 。 这 为 接 下 来 研究 Fe 基 合 金 核反应 
堆 材 料 中 空位 的 形成 位 置 ， 提 供 了 理论 依据 。 


在 不 同位 置 (S1, S2, S3, S4, 01, 02, 03, 04, T1, T2, T3, T4) 


八 面 体 间隙 位 ， 在 


先 倾向 位 于 远离 界 
间隙 位 ， 直 至 靠近 


八 面体 间 际 位 也 填 


满 ， 过 多 的 He 原子 就 会 开始 蔡 代 原 有 的 Cr 原子 与 Fe 原子 ， 取代 他 们 的 位 置 ， 形 成 氨 泡 
或 氨 空 位 复合 体 ， 从 而 产生 微观 缺陷 。 对 于 被 蔡 代 的 Fe 原子 和 He 原子 ， 由 于 脱离 了 原 


有 的 品格 位 置 ， 会 因为 较 大 的 原子 半径 ， 对 原 有 品格 参数 产生 较 大 影 


微观 缺陷 ， 比 如 空位 、 空 位 团 、 团 簇 、 缺 陷 对 、 位 错 、 位 错 环 等 ， 进 


理化 学 性 质 。 
4. 体积 变化 
为 了 更 直观 的 得 昌 


H He 在 Fe-Cr 界面 不 同位 置 的 影响 ， 第 一 性 原理 


啊 ， 在 材料 中 形成 
一 步 影响 材料 的 物 


计算 得 出 He 和 空 


位 在 Fe-Cr 界面 不 同位 置 处 的 体积 如 图 3 所 示 。 
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3 He 和 空位 在 Fe-Cr 合金 界面 不 同位 置 处 的 体积 


Fig. 3 Volume of He and Vacancies at Different Positions of Fe-Cr Alloy Interface 


显而易见 ， 相 比 于 退火 后 没有 微观 缺陷 的 Fe-Cr 合金 界面 , 无 论 是 He 在 替代 位 还 是 


空位 在 蔡 代 位 ， 都 会 使 得 Fe-Cr 界面 的 体积 相对 减 小 ， 且 单 空 位 的 出 现 ， 更 会 明显 降低 
Fe-Cr 界面 的 体积 。 当 He 原子 出 现在 间隙 位 时 ， 会 明显 的 增加 Fe-Cr 界面 的 体积 。 值 得 
注意 的 是 ， 当 He 原子 出 现在 Fe-Cr 界面 的 不 同位 置 时 (Ol, 02, 03, 04, T1, T2, T3, 
T4), T3 位 置 的 体积 比 其 他 位 置 的 体积 都 大 ， 且 体积 的 变化 近乎 是 别 的 位 置 的 两 倍 。 这 

计算 结果 表明 ， 当 He 原子 位 于 靠近 界面 的 Fe 八 面体 间 隐 位 时 ， 会 引起 整个 界面 出 现 
较 大 的 蝇 格 畸 变 ， 这 种 畸变 会 一 定 程度 上 改变 材料 的 形状 ， 从 而 影响 材料 的 物理 化 学 性 


质 。 


针对 于 这 一 位 置 的 特殊 现象 , 我 们 继续 对 He 在 此 位 置 的 结构 模型 进行 了 多 次 精度 更 
高 的 优化 ， 其 优化 结果 都 在 358. 94—359. 2 久之 间 (1A=10 5m) 。 与 此 同时 ， 使 用 p4vasp 
软件 得 出 了 它 的 态 密度 图 以 及 包含 不 同 电子 轨道 的 分 波 态 密度 图 ， 并 与 其 他 八 面体 间 辽 
位 《01, 02, 04) 的 电子 态 密 度 图 进行 了 比较 (图 4)。 
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图 4 He Æ Fe-Cr 合金 界面 不 同 八 面体 间隙 位 01，02，03，04 的 电子 态 密 


坐标 表示 电子 态 密度 


度 图 ， 横 坐标 表示 能 量 区 间 ， 纵 


(03 中 红 绿 蓝 三 条 曲线 分 别 表示 Fe、Cr、He 原子 的 分 波 电子 态 密度 ) 


Fig. 4 Electronic Density of State Charts of He at Different Octahedral Gap Points O1, 02, 03, O4 at Fe-Cr Alloy 


Interface (The three curves of red, green and blue in O3 represent the fractional density of states of Fe, Cr and He 


在 不 同 轨道 


间 
体 


而 


道 与 Fe 原子 和 Cr 原子 的 日 


atoms, respectively ) 


电子 态 密度 能 反应 得 到 电子 在 特定 能 量 区 域内 的 分 布 情况 ， 其 态 密度 图 能 反应 电子 


RZ (Ol, 02, 04) A 


间隙 位 C03) 时 的 波峰 出 现在 能 
置 会 与 周围 的 Cr 原子 和 Fe 原子 有 特殊 的 相互 作用 ， 
Æ Fe, Cr, He 三 种 原子 所 代表 的 红 绿 蓝 三 条 线 都 有 波峰 产生 ， 说 明 He 原子 的 电子 轨 


的 分 布 ， 从 而 反映 出 原子 之 前 的 相互 作用 。 从 图 4 中 可 以 看 出 ， 其 他 八 面体 
Gt -20Ev~-18eV 处 有 一 个 波峰 ， 而 He 在 靠近 界面 的 Fe 八 面 
量 -18eV~-16eV 之 间 。 波 峰 的 偏 移 表示 He 原子 在 此 位 


电子 轨道 出 现 了 杂 化 现象 。 该 峰 是 由 


该 峰 位 不 只 只 是 个 原子 作用 的 结果 ， 


Fe, Cr, He 所 共同 组 成 的 


杂 化 峰 ， 是 Fe、Cr、He 三 原子 之 间 一 种 特殊 的 相互 作用 。 但 是 相 比 于 其 他 八 面体 间隙 位 
来 说 不 稳定 ，He 原子 很 容易 迁移 到 别 的 八 面体 间隙 位 以 及 四 面体 间隙 位 。 


6. 结论 


第 一 性 原理 计算 得 出 ， 当 He 原子 进入 Fe-Cr 合金 界面 时 ， 首 先 倾向 于 向 远离 界面 位 
置 的 Cr 四 面体 间隙 位 置 处 移动 , 当 这 一 位 置 被 完全 占据 , 之 后 依次 是 其 他 四 面体 间隙 位 ， 
八 面 体 间 阶 位 ， 蔡 代位 。 这 一 结论 有 助 于 进一步 研究 核反应 堆 材料 中 的 氮 迁 移行 为 ， 从 
而 研究 材料 的 抗 辐 照 能 力 。 

当 He 原子 进入 靠近 Fe-Cr 合金 界面 的 Fe 八 面体 间隙 位 时 , 会 引起 较 大 的 唱 格 畸变 ， 
这 是 He、Fe、Cr 三 种 原子 在 特定 位 置 相互 杂 化 影响 的 结果 。 这 使 得 界面 的 原子 排列 发 生 
较 大 改变 , 从 而 影响 材料 的 微观 缺陷 结构 , 该 发 现 为 下 一 步 研究 Fe-Cr 界面 内 He 在 Fe 八 
用 体 间 隙 位 时 对 合金 造成 结构 影响 提供 理论 依据 和 研究 思路 。 


T 
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